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ABSTRACT 

The distribution of pairwise, relative peculiar velocities, f(u;r), on small nonlinear 
scales, r, is derived from the Press-Schechter approach. This derivation assumes that 
Press-Schechter clumps are virialized and isothermal. The virialized assumption re- 
quires that the circular velocity, 14 cx M^/^, where M denotes the mass of the clump. 
The isothermal assumption means that the circular velocity is independent of radius. 
Further, it is assumed that the velocity distribution within a clump is Maxwellian, 
that the pairwise relative velocity distribution is isotropic, and that on nonlinear scales 
clump-clump motions are unimportant when calculating the distribution of velocity 
differences. 

Comparison with TV-body simulations shows that, on small nonlinear scales, all 
these assumptions are accurate. For initially scale-invariant Gaussian density fields, the 
pairwise velocity distribution evolves in a self-similar manner. This is consistent with 
other analytic expectations, and with the distribution measured in relevant TV-body 
simulations. For most power-spectra of interest, the resulting line of sight, pairwise, 
relative velocity distribution, f(uj-), is well approximated by an exponential, rather 
than a Gaussian distribution. 

This simple Press-Schechter model is also able to provide a natural explanation for 
the observed, non-Maxwellian shape of f(v), the distribution of peculiar velocities. 

Key words: galaxies: clustering - galaxies: evolution - galaxies: formation - cosmol- 
ogy: theory - dark matter. 



INTRODUCTION 



When quantifying the relation between real and redshift 
space galaxy catalogs, the distribution function of pairwise 
peculiar velocities is of fundamental importance. In the lin- 
ear regime the effects of redshift induced distortions can 
be calculated explicitly (e.g. Peebles 1980). On large, linear 
scales, the correlation function in redshift space is closely re- 
lated to the correlation function in real space (Kaiser 1986; 
Hamilton 1992). In the nonlinear regime, also, the correla- 
tion functions in real and redshift space may be easily re- 
lated, provided that on small, nonlinear scales the clustering 
is stable so that, on average, the pairwise motions (on these 
small scales) exactly cancel the Hubble expansion. However, 
to relate the two correlation functions requires knowledge of 
the distribution of pairwise relative peculiar velocities (Pee- 
bles 1980, §76). It is known from the linear theory that if the 
initial density field is Gaussian, then the pairwise relative ve- 
locity distribution that obtains on large, linear scales must 
also be Gaussian (e.g. Fisher 1995). In contrast, on smaller 
nonlinear scales, observations (Peebles 1976) and A^-body 
simulations (Efstathiou et al. 1988; Fisher et al. 1994) sug- 



gest that the pairwise relative velocity distribution, when 
viewed along the line of sight, is better approximated by 
an exponential, than by a Gaussian distribution. This pa- 
per shows that the accuracy of the exponential distribution 
on nonlinear scales can be easily explained by the Press- 
Schechter theory of nonlinear clustering. 

Section 2 shows how to calculate the distribution of 
pairwise, relative peculiar velocties. It restricts attention to 
the highly nonlinear regime, since when pair separations are 
small the distribution of pairwise velocities should be insen- 
sitive to clump-clump motions. Essentially, this is because 
particles that are separated by small distances are most 
likely to be in the same clump, so the motion of the clump 
as a whole is irrelevant when computing the pairwise rela- 
tive velocity. The calculation assumes that Press-Schechter 
clumps are virialized and isothermal, and that the velocity 
distribution of particles within each clump is Maxwellian. 
For initially scale-invariant power-spectra, the distribution 
of pairwise relative peculiar velocities that results evolves 
self-similarly, in accordance with the prescription discussed 
by Davis fe Peebles (1977), and tested by Ffstathiou et al. 
(1988). For all power-spectra of current interest, the distri- 
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bution of line of sight, pairwise, peculiar velocities is well fit 
by an exponential, rather than by a Gaussian distribution. 

Section 3 shows how to use the model developed in Sec- 
tion 2 to calculate the peculiar velocity distribution function 
(as opposed to the distribution of velocity differences) using 
the Press-Schechter approach. This distribution is the ana- 
logue of the Maxwell-Boltzmann distribution that describes 
an ideal gas. Unlike the pairwise, relative peculiar velocity 
distribution computed in the previous section, this distri- 
bution can only be derived assuming some model for the 
motions of clumps towards each other. Section 3 presents 
a simple, ad hoc model for this clump-clump velocity dis- 
tribution. In this respect, the results are Section 3 are less 
general than those of Section 2, since they dependend on as- 
sumptions about the motions of clumps towards each other. 
The peculiar velocity distribution function that results is 
well fit by a functional form, used previously by Saslaw et 
al. (1990), that fits A^-body simulations of clustering from 
Poisson initial conditions well. It is interesting that this dis- 
tribution also appears to be consistent with the measured 
peculiar motions of galaxies (Raychaudhury fe Saslaw 1995) 
that are available at present. 

Section 4 summarizes the results and discusses the re- 
lation between the assumption that clustering on nonlin- 
ear scales is statistically stable (i.e., that, on average, the 
pairwise velocities cancel the Hubble flow exactly), and the 
Press-Schechter model that is developed and extended in 
this paper. 



2 THE DISTRIBUTION OF PAIRWISE, 
RELATIVE PECULIAR VELOCITIES 

This section develops a model for the distribution of pair- 
wise, relative peculiar velocities for pairs that are separated 
by small distances. 



2.1 The model 

Let f{v) dv denote the probability that the peculiar velocity 
of a particle is in the range dv about v. Also, write the 
peculiar velocity of the ith particle as the sum of two terms: 



Up, + Vc 



(1) 



where the first term on the right describes the motion of the 
particle with respect to the center of mass of the clump 
of which it is a member, and the second term describes 
the motion of the clump as a whole. Then the pairwise 
velocity distribution for pairs separated by the vector r, 
f{u;r) = f{vi — V2',r), requires knowledge of the joint ve- 
locity distribution, f{vi,V2), where particles 1 and 2 are 
separated by r. This with equation (1) shows that, in gen- 
eral, to compute f{u;r ) requires knowledge of the pairwise 
velocity distribution of clumps. Usually, these clump-clump 
velocities are not known. However, if we restrict attention 
to separations that are sufficiently small that both mem- 
bers of most pairs are likely to be in the same clump, then 
f{u;r) can be obtained without any knowledge of the mo- 
tions of clumps, since then Uij = Dp, — Vp^ for most pairs. 
So, provided we restrict attention to separations that are 
small compared to the dimensions of a typical clump, the 
distribution of pairwise velocities can be obtained without 



any knowledge of the motions of clumps. Clearly, however, 
the distribution of pairwise velocities, even on these small 
scales, depends on the distribution of velocities within a 
given clump, and on the distribution of clump sizes. So, for 
small separations, the pairwise velocity distribution can be 
written as 



f{u;r) 



p{u\M;r)p{M;r)dM, 



(2) 



where p{u\M; r) is the probability that the pairwise velocity 
difference for particles separated by r in a clump of mass M 
is u, the second term gives the probability that the pair 
is within a clump of mass M, and the integral is over all 
clumps that are larger than \r\. In particular. 



p{M;r)dM 



iVpa.rs(M;r-)dM 
/ 7Vpa.rs(M;r)dM' 



(3) 



where the denominator denotes the total number of pairs 
that have the separation vector r. For small values of \r\, 
this total number of pairs may be approximated by assum- 
ing that both members of each pair are in the same clump, 
so that the integral is over all clumps that are larger than 
\r\. In what follows, we will assume that the distribution of 
u = Vi — V2 is isotropic. On the small scales of interest, this 
is likely to be a good approximation (e.g. Efstathiou et al. 
1988). This means that p{u\M; r) and p{M;r) may be re- 
placed with p{u\M;r) and p{M;r), respectively. With this 
assumption of isotropy, f{u; r) will be computed by using 
the Press-Schechter description of nonlinear gravitational 
clustering (Press fe Schechter 1974; Bond et al. 1991) to 
construct models for both p{u\M; r) and p{M; r). 

In the Press-Schechter approach, the number density of 
clumps of a given mass can be computed directly from the 
statistics of the initial density field. These Press-Schechter 
clumps are assumed to virialize when they are ~ 178 times 
more dense than the background density. This fact will be 
useful when computing the number of pairs as a function of 
separation, p{M; r). A virialized clump of mass M may be 
assigned a circular velocity, Vc, in accordance with the pre- 
scription suggested by Narayan fe White (1988), and used 
later by White fe Frenk (1991). Namely, the assumption that 
Vc = GM I fvir , where Vc defines a circular velocity that can 
be associated with the clump of mass M , and rvir denotes 
the radius of the clump (i.e., pvir = 178/9b = 3Af/47rrJ;j.), 
provides a relation between the circular velocity Vc and the 
mass: Vc (x M^^^. Assuming that these virialized clumps are 
also isothermal means that Vc is approximately independent 
of position within the clump. A^-body simulations show that 
the Vc oc M^l^ scaling is quite accurate (e.g.. Bond fe Myers, 
in preparation; but see discussion at the end of this paper). 
The simulations of Efstathiou et al. (1988) show that, within 
a given clump, Vc is, indeed, approximately independent of 
radius. 

For an isothermal sphere Vc is two-thirds of the three 
dimensional velocity dispersion (cf. Section 4.4.3(b) in Bin- 
ney fe Tremaine 1987; Appendix B in Lacey fe Cole 1993). 
Therefore, in this model, particles in a clump of mass M are 
assigned velocities from a Maxwellian distribution that has a 
three-dimensional mean square velocity given by ?>Vc{M)l2. 
The choice of a Maxwellian for the velocity distribution 
within each clump is rather ad hoc; although it is suggested 
by the virial assumption, it is not required. It may be that 
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one of the self-consistent velocity distributions derived by 
Evans (1994) for density distributions that vary as power- 
laws of the radius, is more appropriate. The Maxwellian is 
chosen here simply to provide some insight into the shape of 
the velocity distribution when it is averaged over all Press- 
Schechter clumps. 

Thus, in this model, the distribution of pairwise veloc- 
ities on small scales is given by calculating the distribution 
of the vector u = Vi — V2, where Vi and V2 are Maxwellian 
vectors. For small separations, it is very likely that both 
members of a given pair are members of the same clump. 
So, the assumption that the clump is isothermal means 
that the Maxwellian distribution that describes Vi also de- 
scribes the distribution of 1)2. Therefore, the distribution of 
u = Vi — V2 is also Maxwellian, since calculating the dis- 
tribution of Vi — 1)2 reduces to computing the distribution 
of {vxi — Vx2,Vyi — Vy2,Vzi — Vz2), whcrc the VijS are all 
Gaussian distributed. Since it is obtained from the differ- 
ence of two Gaussian random variables that have the same 
mean and variance as each other, the mean square of this 
pairwise Maxwellian velocity distribution is ?>V^ [M] (this 
is essentially Problem 7-3 on p. 485 in Binney fe Tremaine 
1987). This follows because the three dimensional disper- 
sion of the Maxwellian that describes the motions of the 
particles themselves is ?>V^{M)l2, and we are noting explic- 
itly that the pair belongs to a clump of mass M . Thus, the 
assumption that the Press-Schechter clumps are virialized 
isothermal spheres implies that the distribution p(u\M; r) of 
equation (2) is Maxwellian, with a dispersion that is related 
to the mass of the clump. 

The distribution of pairs, p{M; r), can be obtained from 
the assumption that all Press-Schechter clumps are trun- 
cated, singular isothermal spheres having the same density 
(~178 times the background density). The Appendix shows 
that if the distribution of particles is isotropic, then the num- 
ber of pairs separated by a distance r that are in a clump of 
mass M and radius R scales as 7Vpairs(M;r) oc M'^/R^, at 
least to first order. Since all clumps have the same density, 
oc M, so that 7Vpairs(M; r) oc M^'\ Thus, 



p{M;r)dM 



M^l^ n(M) AM 
Im ^^'^ n(M)AM' 



(4) 



where n(M) is the number density of clumps of mass M. 
Note that, for consistency, the integral in the denominator 
is only over those clumps whose diameters, 2R, are larger 
than the pair separation r. Since all clumps have the same 
density, this implies that Afmin c>c -Rj^in, with -Rmin = rl2. 
Strictly speaking, equation (4) for the pair distribution is 
correct only for pair separations that are small compared to 
the clump radius. The exact expression, given in the Ap- 
pendix, is complicated. For the sake of simplicity (and be- 
cause it will turn out that the results are not very sensitive 
to this approximation), we will continue to use the simpler 
equation (4) in the remainder of this section. 

Since n(M) is given by the Press-Schechter distribu- 
tion, f{u; r) can now be obtained by computing the inte- 
gral in equation (2). The assumption of isotropy allows us 
to obtain the line of sight distribution from equation (2). 
Since the velocity dependence is contained entirely in the 
Maxwellian p(u\M;r) term, it is sufficient to integrate this 
Maxwellian over all directions perpendicular to the line of 



sight. The Maxwellian reduces to a Gaussian, so that the 
line of sight distribution is a Gaussian (with dispersion 
cr^ = 2(73£, = ZV^ oc 3M^''^) convolved with the p(M;r) 
distribution of equation (4): 



/(«; '•) 



p{M;r)AM, 



(5) 



where u now denotes the line of sight, rather than the three 
dimensional velocity. Once n(M) in p{M; r) has been spec- 
ified, f{u;r) may be obtained by integrating equation (5). 
Equation (5) is the main, general result of this paper. It 
shows that the line of sight pairwise velocity distribution is 
essentially a weighted sum over Gaussians having a range of 
dispersions, where the weighting factor for each Gaussian is 
related to the Press-Schechter multiplicity function, and to 
the distribution of particles within a Press-Schechter clump. 



2.2 Relation to BBGKY scaling solution 

This subsection considers a useful, specific example of the 
model derived above; the evolution of the pairwise velocity 
distribution as an initially scale free density field clusters 
gravitationally. 

For an initially scale free field, Davis fe Peebles (1977) 
argue that the BBGKY equations can be transformed into 
a series of equations for the moments of the pairwise ve- 
locity distribution. They show that the A^-point correlation 
functions admit a similarity transformation of the form 



s = where 



3(3 -I- n) 



(6) 



where n is the slope of the power-spectrum of the initially 
scale invariant Gaussian density field [i.e., the initial power 
spectrum is P{k) oc k"], t denotes cosmic time, a denotes 
the expansion factor at that time (a oc t^^^ in a Universe 
that has critical density), and a; is a comoving distance. 
They also argue that if equation (6) is satisfied, then the 
BBGKY hierarchy requires that the relative peculiar veloc- 
ities be rescaled by the relation 

(7) 



(e.g. Efstathiou et al. 1988). We can compare this scaling 
solution to the model described in the previous subsection. 

Recall that GM/r,,, = = GM/(3M/47rAnl/9b)^''^ 
where Ani ~ 178 and ph is the background density, 

and that = where cr^ is the three dimensional 

velocity dispersion of the pairwise velocity distribution. In 
a Universe with critical density (i.e. 3H^ = SwGph), these 
relations imply that 



M 



G 



a/2 



47rGAni/9b J G S^/^GHAIC" 



(8) 



Equation (8) defines a, and a oc a^^^ since pb oc a ^ . This, 
in equation (5), gives 



/(«; '•) 



where 



(9) 
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„(M)dM=^ f^)"*'e-(^)"V^) (10) 

is the Press-Schechter distribution for an initially scale free 
field, 



(11) 



and Cn is a constant that depends on the slope n and am- 
plitude of the initial power spectrum (e.g. Efstathiou et al. 
1988). Now define 



u = (a/Mt)^^^ u (X. 



(2/(3 + n) ^4/3(3 + n) 



(12) 



where A was defined in equation (6). So u scales similarly 
to the Davis fe Peebles solution. Also, set = M/Mt. Then 



/(«; f) 



/(«;»•)— = f{n;r) 

du \ a ' 



J » tt + 3 



t + 3 

-1/ 3 / n + 3 \ dv 
V 3 J 1/ 



7 1-+ — f 

M 2 ^ n + 3 ' mm 



,(13) 



where the denominator in the final expression is the com- 
plimentary incomplete Gamma function. This shows ex- 
plicitly that if J'min is constant, then f{u; f ) evolves self- 
similarly. However, recall that Afmin is related to -Rmin and 
f/2, where r is the separation scale. This means that 



i^min 



M"min _ 47rAnl/qo I -Rmin/a 
M. ~ 3Cn U*/3(3 + n) 



(14) 



Therefore, requiring that J'min remain constant implies 
that _Rmin/a, the comoving cutoff distance, must scale as 
^4/3(3+n) Comparison with equation (6) shows that this is 
just what is required by the Davis-Peebles solution. Thus, 
this model for the pairwise velocity distribution satisfies the 
Davis-Peebles scaling solution of the BBGKY hierarchy. 
Having shown that this Press-Schechter model for the 



pairwise velocity distribution on small scales has the correct 
scaling, we now consider the shape of this distribution. First, 
consider the limit of very small pair separations (so that 
Afmin 0). The pairwise velocity distribution becomes 



/(«) 



f CK 

Jo 



e '■' ■'2/3 J/ 6 e 



poo "> 

Jo " 



■ ( n+3 \ dv 

\ !. I V 

+ dx^ 

\ 2 ) 



(15) 



phere x 



,2/3 



r (i -h 

V2 ^ n + 3 / 

{MIM^fl''. In this limit, the pairwise 



dispersion is 

1 r (i -2- 

J- V2 ^ n + 3 



V2 ^ n+3/ 

Equation (16) shows that for n 
rms line of sight velocities are a, 
and ^1/3, respectively. 

When n = —1, then equation (15) can be calculated 
analytically: 



(16) 

= 1,0, and —1, the scaled 
'{m2) = VO.2466, VO.2617, 



/(«) 



.3 «2 _^ dx/ 
X e — — 



exp 



(17) 



The final expression follows from equation (3.325) in Grad- 
shteyn fe Ryzhik (1994). Equation (17) shows that when the 
initial fluctuation spectrum is scale free with slope n = —1, 
then in the limit of vanishingly small pair separations, the 
probability that the pairwise velocity is in the range du 
about u is given by an exponential distribution. 

Moreover, consider those initial conditions that have 
Press-Schechter n{M) distributions that are not vastly dif- 
ferent from the n = —1 shape. In the limit of very small pair 
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separations, the exponential will provide a good approxima- 
tion to f{u) for these other initial conditions also. To illus- 
trate this, Fig. 1 shows f{u) distributions for ra = 1,0 and 
— 1, in the limit where J'min = 0. The solid line in each panel 
shows equation (15). The corresponding value of the mean 
square pairwise velocity (eq. 16) is also shown. Dot-dashed 
and dotted lines show exponential and Gaussian distribu- 
tions that have the same dispersion as the solid curves. In 
all three panels, the solid curves (eq. 15) are well approxi- 
mated by the dot-dashed curves (exponential distributions). 

When the cutoff, J'min, is not vanishingly small, then 
the shape of f{u; f) will depart from the solid curves shown 
in the figure. The extent of these departures may be esti- 
mated as follows. A given value of J'min = Afmin/Af» means 
that clumps less massive than Mrain/Mt = (-Rmin/'"*)^ do 
not contribute to f{u; f). Now, smaller clumps have smaller 
circular velocities, so they contribute to f{u; f ) primarily 
when |m| is small. Therefore the existence of a nonvanishing 
J'min has two related effects. The first is that, for a given 
I'min, the departures from the J'min = shape will be maxi- 
mal near the small |m| peak of the distribution. This implies 
the second effect: As i'min increases, the dispersion will also 
increase. 

However, this model for f{u; f) (eq. 9) only applies 
for small values of i'min. To see this, recall that r, is the 
scale on which the mass variance in spherical volumes is 
X Sl/2 X (1.68/a)V2 (e.g., Efstathiou et al. 1988). So, 
r, should be slightly larger than ro, the correlation length, 
which is approximately the average Press-Schechter clump 
size (Peebles 1980). Since this model is only accurate on 
scales that are smaller than the average clump size, i'min = 
(-Rmin/j"*)^ ~ (-Rmin/''o)^ should always be less than, say, 
(1/2)^ ~ 0.1. For these small values of i'min, the tails of 
f{u; f ) should remain unchanged from their i'min = val- 
ues. To see the effects of i'min directly, equation (9) must be 
integrated numerically. Fig. 2 in the next subsection shows 
the result of integrating equation (9) for a range of values of 
ra, and for representative values of the cutoff. 



2.3 Comparison with A^-body simulations 

The Press-Schechter model for f{u; f ) developed in the 
previous subsection may be compared with that measured 
in the A^-body simulations of gravitational clustering from 
scale free initial conditions studied by Efstathiou et al. 
(1988). Efstathiou et al. define scaled distances and veloci- 
ties using the transformations given in equations (6) and (7). 
They define so as that value of the similarity parameter, s, 
at which the correlation function is unity, and present most 
of their velocity distribution results in terms of the variable 
mefwd/so, scaled to a = 1. They find that the scaled veloc- 
ity distributions that they measure in their simulations have 
roughly exponential tails for the values of ra they consider. 
However, the peaks of the distributions are more rounded 
than an exponential (their Fig. 6). 

Since the model developed in the previous subsection 
also produces exponential distributions (Fig. 1), we must 
check that it also gives the measured dispersions. To do so, 
we must relate the scaled velocity, u, defined in equation (12) 
to mefwd , the scaled variable used by Efstathiou et al. Using 
equations (8), (11), (12), and the fact that ao = 1 and Cn = 



0.8,0.71, and 0.53 for ra = 1,0, and —1, respectively (their 
eq. 15), yields 



a 

m7 



1 / 2 \ 6 / 87r 

3 [a^J V3a 



27r 2 at-*-! 



Since Ani = 178, this means that 



mefwd 0.236 C, 



1/3 



So 



Soto 



,(18) 



(19) 



Their Fig. 4 shows that sot^/L X 0.048/27r, 0.054/2 TT and 
0.065/27r, for ra = 1,0, and —1, respectively. These values, 
with their Fig. 6, suggest that to = 2/3. In the limit of very 
small pair separations we can use equation (16) to estimate 
the variance of the model distributions: 



)> 0.236 C, 



So 




(20) 



Equation (20) implies that the scaled one dimensional rms 
pairwise dispersion in the simulations should take the val- 
ues 1^ {y?^ I So = 2.2,2.0, and 1.7, for ra = 1,0, and —1, re- 
spectively. These values differ slightly from those measured 
in the simulations, where the corresponding values are 2, 
1.8, and 1.6 (the small s limit in Fig. 5 of Efstathiou et al. 
1988). Some of this discrepancy may be due to the fact that 
in the simulations, virialized clumps collapse by a factor of 
1.8 rather than the virial factor two (Hamilton et al. 1991, 
who suggest that some of this may be due to force softening 



effects), so that A^j'^ should be replaced by 0.9A^j'''. How- 
ever, some of the discrepancy may also be due to the fact 
that not all clumps in the simulations are truncated singu- 
lar isothermal spheres. Whatever the reason, these values 
for the dispersion are close enough to the A^-body results to 
consider the model further. 

Having shown that the Press-Schechter model has the 
correct scaling and approximately the right pairwise veloc- 
ity dispersion, we now examine the shape of the distribu- 
tion for nonzero values of the cutoff. Recall that this shape 
is given by equation (13) which must be integrated nu- 
merically. For ease of comparison with the A^-body simu- 
lations, a representative value of the cutoff, i'min, is chosen 
as follows. In their Fig. 4, Efstathiou et al. show the cor- 
relation functions measured in their simulations, and indi- 
cate the scales xioo and xo at which { = 100 and { = 1. 
This allows us to estimate xio- Their figure suggests that 
logj^o (a^o/a^io) ~ 0.3,0.38 and 0.45 for ra = 1,0, and — 1, 
respectively. The last paragraph of the previous subsec- 
tion argued that i'min ~ (-Rmin/''o)^ ~ {x/xoY- So, setting 
I'min ~ {xio/xoY givcs au indication of the shape of the 
velocity distribution for particles separated by the scale cor- 
responding to { ~ 10. 

Fig. 2 shows /(mefwd/so) as a function of mefwd/^h, 
when ra = 1, 0, and —1 for this choice of i'min. Velocities are 
in units of the Hubble velocity across the Efstathiou et al. 
(1988) computational volume, and have been rescaled using 
equation (19) to a = 1, for ease of comparison with their 
Fig. 6. (On the scale corresponding to s, the scaled value 
of the Hubble velocity is = 2s/3.) Each panel shows 
the value of ra and the corresponding rescaled value of the 
mean square pairwise velocity. Solid lines show the distri- 
bution that is obtained by integrating equation (13) with 
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u/vrr 'i^/'^ff u/v 



Figure 2. Examples of the distribution function of pairwise relative peculiar velocities on scales where ^ % 10. Velocities are in units 
that simplify comparison with the scaled /(«efwd /^o) distribution measured in A^-body simulations (Fig. 6 of Efstathiou et al. 1988). 
In each panel, solid line shows equation (13), scaled to the A^-body variables as described in the text (eq. 19). The corresponding scaled 
value of the mean square pairwise velocity is also shown. Dotted lines show Gaussian distributions that have the same dispersions as 
the solid curves. Dot-dashed lines show exponential distributions with dispersions given by equation (20); for these values of n, these 
dispersions are slightly less than those of the solid curves. 



= 0.126,0.075, and 0.045 for ra = 1,0, and -1, respec- 
tively. Dot-dashed lines show exponential distributions with 
the same dispersion as the solid curves, had the solid curves 
been computed with J'min = 0. Comparison of the shapes of 
these curves with those in Fig. 1 shows that the effect of the 
cutoff is most pronounced near the peak of the distribution; 
the tails are hardly affected. This is consistent with the dis- 
cussion at the end of the previous subsection. Dotted lines 
show Gaussian distributions (corresponding to Maxwellians 
in three-dimensions) with the same dispersions as the solid 
curves. 

To conclude this subsection, we note that the solid curve 
in the leftmost panel of Fig. 2 (ra = 1, on scale correspond- 
ing to { Si 10) resembles the leftmost panel of Fig. 6 of 
Efstathiou et al. closely. Some of the small differences be- 
tween the two figures may be due to the fact that whereas 
Efstathiou et al. show the ra = 1 distribution for { Si 10, 
the scale for Fig. 2 is only approximately where { Si 10. 
Moreover, recall that Fig. 2 was constructed using an ap- 
proximate expression for p{M; r) (cf. eq. 4). Using the exact 
expression (given in the Appendix) makes the two figures 
even more alike. In any case, the figures are similar enough 
that we conclude that this Press-Schechter model for the 
pairwise velocity distribution is in good agreement with the 
A^-body simulations. 

2.4 The effects of discreteness 

There is one complication when comparing these Press- 
Schechter models of clustering from initially Gaussian den- 
sity fields with the A^-body simulations. It arises from the 
fact that, in the Press-Schechter description, the probability 
of forming arbitrarily small clumps is not negligible. Since 
the velocity dispersion within a clump is related to its mass. 



this means that the probability of having arbitrarily small 
velocities is also not negligible. The A^-body simulations, on 
the other hand, only have a limited dynamic range; they 
cannot resolve clumps that are smaller than some minimum 
size. In terms of the pairwise velocity distributions, this 
means that the simulations are only sensitive to clumps that 
have velocity dispersions greater than some minimum value. 
When comparing this Press-Schechter model with simula- 
tions, it is not clear what effect this discreteness will have 
when choosing the lower bound on the integration variable 
in equation (13). 

The importance of discreteness in the Efstathiou et al. 
(1988) simulations can be estimated as follows. Assume that 
the particles in the simulations constitute a distribution of 
points that has the same correlation function as the con- 
tinuous density field that the particles are being used to 
simulate. One way to do this is to assume that the parti- 
cles in an A^-body simulation represent that point distri- 
bution which is obtained when the Poisson sampling pro- 
cess discussed by Layzer (1956) and by Peebles (1980) is 
applied to the continuous density field. Indeed, the assump- 
tion of just such a relation between the discrete particles in 
the simulations and the underlying continuous density field 
is necessary to justify comparing statistical quantities such 
as the discrete particle-particle correlation function in these 
simulations with theoretical calculations of the evolution of 
the two-point corrrelation function of the continuous density 
field. 

If the Poisson sampling process is an accurate model of 
the role of discreteness in the simulations, then the way in 
which the Press-Schechter mass functions are modified by 
the discreteness is calculable. Although, by assumption, all 
clumps have exactly the same overdensity (~178 times the 
background density), the calculation is made easier by as- 
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Figure 3. Multiplicity functions that are obtained after applying the Poisson sampling process of equation (21) to the Press— Schechter 
n[M) distributions, for three choices of n, and a range of epochs. Plots show the fraction of particles at each epoch that belong to groups 
with Nm members, where 2^~^ < < 2"^. Solid lines show the discrete clump multiplicity functions (obtained from eq. 21) for a 

range of expansion factors (all but the last one as shown in Fig. 8 of Efstathiou et al. 1988). Dotted lines show the corresponding curves 
for the continuous distribution M'^n[M) logg 2. 



suming that, to first order, all clumps have the same volume. 
For the initially scale free distributions studied by Efstathiou 
et al. (1988), this means that, to a good approximation, 
the Poisson sampled discrete Press-Schechter functions are 
given by 



V{N) 



n{M) dM, 



(21) 



where riiN) is the probability that a Press-Schechter clump 
has N associated particles, so that the denominator on the 
left is the average clump size. When ra = this integral is 
analytic: 



V{N) 



1 



E^viN) r(7v + i)r(i) VI + M. VI + M. 

where the sum is over all > 0. Thus, 



— (a/ttm:-!) 



(22) 



(23) 



En=i Nv{N) M.. 

which has been simplified using the definition (2ra)!! = 2" ra!, 
the identity r(ra + |)/r(|) = (2ra - l)!!/2", and the Tay- 
lor series expansion of (1 - xfl'^. Requiring that ii{N) be 
normalized to unity implies that 



1 r (TV - f 



2 r (TV + 1) r (i) V 1 + M. ; ^i + m, - i 



(24) 



For other values of ra, equation (21) must be integrated nu- 
merically. 

The panels in Fig. 3 show the functions that are ob- 
tained after applying the Poisson sampling process of equa- 
tion (21) to the Press-Schechter ra(Af) distributions, for a 
range of choices of ra. The solid lines in each panel show the 
resulting discrete clump multiplicity functions (eq. 21) for 
a range of expansion factors (the first six of the last seven 



expansion times in the Efstathiou et al. simulations), and 
are chosen for comparison with Fig. 8 of Efstathiou et al. 
(1988). (As Mt increases, the curves peak at larger values of 
m.) For comparison, the dotted lines show Af^ra(Af) log,, 2, 
where ra(Af) is the continuous Press-Schechter distribution, 
at these same epochs. 

For each ra, the solid (discrete) and dotted (continu- 
ous) distributions have similar large clump tails, but dif- 
fer increasingly as the clump size decreases. This difference 
decreases as the characteristic mass Af, increases. This be- 
haviour is easily understood since, as the simulation evolves, 
the characteristic mass increases, so the mass of a single par- 
ticle as a fraction of the characteristic mass decreases. Thus, 
the effects of this lower mass cutoff, this discreteness, be- 
come less pronounced as the simulation evolves. The multi- 
plicity functions in Fig. 3 are similar to those shown in Fig. 8 
in Efstathiou et al. (1988). (Applying the Poisson sampling 
process to the modified mass functions given by Jedamzik 
1995 increases the similarity between Fig. 3 here and the 
Efstathiou et al. figure significantly.) Thus, Fig. 3 suggests 
that the Poisson sampling process does, indeed, provide a 
good description of the relation between the particles in the 
A^-body simulations and the underlying idealized continuous 
density field. 

We have established that the Poisson sampling process 
of equation (21) provides a self-consistent way of modelling 
discreteness effects on the evolution of structure in the A^- 
body simulations. So, we can use this sampling model to 
estimate the importance of discreteness on the distribu- 
tion of pairwise velocities. This is done by summing over 
the discrete r](N) distribution, rather than integrating over 
the continuous ra(Af) distribution, when computing the ana- 
logue of equation (9). In the limit of vanishing separation, 
the velocity distributions obtained using the r](N) distribu- 
tions resemble those in Fig. 1 closely, with the agreement 
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increasing at later times. This is expected since the continu- 
ous and the discrete multiplicity functions are quite different 
at early times (Fig. 3). For example, when ra = 0, then the 
fractional difference between the mean square velocity in 
the continuous and the discrete cases is less than one per 
cent for expansion factors a > 4. The curves obtained using 
viN) have slightly more rounded peaks than those in Fig. f . 
Since the discreteness provides a lower bound to the mass, 
this slight rounding of the peaks is consistent with the dis- 
cussion at the end of Section 2.2. We conclude, therefore, 
that the effects of discreteness on the shapes of the pairwise 
velocity distributions in these simulations are not significant. 

Measurements of /(m; r) from observations of galaxy pe- 
culiar velocities (Peebles f976) and A^-body simulations of 
gravitational clustering (Fisher et al. f994) show that, at 
least on highly nonlinear scales, an exponential distribution 
is a good approximation to the exact f{u;r) distribution. 
This section suggests that the Press-Schechter approach is 
able to provide a simple, natural explanation for this agree- 
ment. 



that relates the dispersion of s to the dispersion of c is inde- 
pendent of N. Finally, set the constant factor by requiring 
that the dispersion of d = s + c, when averaged over all 
particles in all clumps, gives the three-dimensional velocity 
dispersion, cr^. 

As was the case in Section 2, the resulting distribution 
(now simply of peculiar velocities, v = s + c, not pairwise 
peculiar velocities, u) is given by convolving a Maxwellian 
with the distribution of circular velocities that is prescribed 
by the Press-Schechter approach. The resulting three di- 
mensional velocity distribution is given by a relation like 
equation (5), with u replaced by v, and where the Press- 
Schechter mass function is weighted by the number of parti- 
cles, rather than by the number of pairs. That is, the proba- 
bility f{v) that a particle has velocity in the range dv about 
V is 

f CO 

f{v)= / p{v\M)Mn{M)dM, (25) 
Jo 

where p{v\M) is a Maxwellian with dispersion related to M, 
and n{M) is the distribution of clump masses. 



3 THE DISTRIBUTION OF PECULIAR 
VELOCITIES 

This section considers the distribution f{ v) of peculiar veloc- 
ities, which is the analogue of the Maxwell-Boltzmann veloc- 
ity distribution that obtains in a perfect gas. It is, of course, 
different from the pairwise velocity distribution, f{u), con- 
sidered in the previous section. 

3.1 The model 

To begin, we summarize the relevant features of the model 
developed in the previous section. Press-Schechter clumps 
are assumed to be virialized singular isothermal spheres. 
The virialized assumption requires that the circular velocity, 
Vc (X M^^^ . The isothermal assumption means that the cir- 
cular velocity Vc is independent of radius. Within a clump, 
the distribution of velocities relative to the center of mass 
is assumed to be Maxwellian. Whereas in the previous sec- 
tion clump-clump motions were irrelevant, they are the main 
source of uncertainty here. To compute f{v) we need a model 
for these clump motions. 

Consider the following simple model for these clump 
motions. Although it is not fully self-consistent, and a more 
detailed treatment, such as that developed by Bond fe My- 
ers (in preparation) may be more appropriate, this simple 
model will serve to provide some insight into the shape of the 
resulting velocity distribution. First, consider those clumps 
with exactly N identical member particles. Let cj denote 
the velocity of the j'th such clump, i.e., the motion of the 
center of mass of the j'th clump. Let Sij denote the velocity, 
relative to the center of mass of the clump, of the ith particle 
in the j'th clump. As in the previous section, assume that 
the distribution of s is Maxwellian (with a dispersion that is 
related to N). In addition, assume that the distribution of c 
is also Maxwellian, and that its dispersion is related to the 
dispersion of s by some constant factor. Since s and c are 
both Maxwellian, 5 -|- c is also Maxwellian (e.g. Binney fe 
Tremaine f987, p. 485). Now allow for the Press-Schechter 
distribution of clump sizes N, and assume that the constant 



3.2 Comparison with A^-body simulations 

The distribution of peculiar velocities has not been stud- 
ied as extensively as has the distribution of pairwise veloc- 
ities. However, an extensive study, using A^-body simula- 
tions, of clustering from an initially Poisson distribution has 
been made by Itoh, Inagaki fe Saslaw (f988; f993), and by 
Saslaw et al. (f990). In addition to measuring various prop- 
erties of spatial clustering, they also measure the velocity 
distribution. Therefore, we will compare the model (eq. 25) 
with the peculiar velocity distribution they measure in their 
simulations. To do so, we must compute the Poisson Press- 
Schechter mass multiplicity function. 

The Press-Schechter description of clustering from an 
initially Poisson distribution has been derived (Epstein f 983; 
Sheth f 995). The probability that a Poisson Press-Schechter 
has A^-particles is given by the Borel distribution: 

and 8c ~ f.68/a is the usual Press-Schechter critical den- 
sity that is necessary for collapse at the epoch a. This dis- 
tribution of clump masses is in good agreement with the 
spatial distribution of particles in the Itoh et al. A^-body 
simulations (Sheth f995). At any epoch, the Poisson Press- 
Schechter model suggests that the spatial distribution of the 
particles is described by 

P(N,nV) = ^^^j^ [7V(f-6)+A^6]^"'e-^(^-^)-^\(27) 

where P{N, nV) denotes the probability that a randomly 
placed cell of size V contains exactly A^ particles, n is the 
average number density of particles, A^ = nV denotes the 
average number of particles in a y-sized cell, and b increases 
monotonically from zero (initially) to some larger value (not 
greater than unity) as the clustering develops (Sheth f995). 
Itoh et al. confirm the accuracy of equation (27) for the spa- 
tial counts in cells distribution. (Incidentally, Sheth, Mo fe 
Saslaw f 994, and references therein, show that eq. 27 also de- 
scribes the spatial distribution of optical and infrared galax- 
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Figure 4. Example of the distribution function that corresponds to the Poisson Press— Schechter model. Histogram shows equation (25) 
when it is required to have the same velocity dispersion as the simulations when 10^ particles have collapsed into Press— Schechter clumps, 
in such a way that the spatial distribution of the particles is described by equation (27) with b = 0.75. Dotted line shows a Maxwellian 
that has the same dispersion as the histogram. Dashed line shows the best fit of equation (28) to the A^-body histogram; top right 
corner shows the values of the parameters, a,/3, and B, that provide the best-fit to the histogram shown. These best-fit parameters are 
very similar to those for the A^-body curve. Velocities are in "natural units" for ease of comparison with the Itoh et al. (1993) A^-body 
simulations. This is why /3, rather than a^, is shown. 



ies well.) Given that the spatial distribution of the particles 
is well described by the Press-Schechter approach, we can 
ask if the velocity distribution is also. 

Fig. 4 shows a comparison of the Press-Schechter model 
for f{v) with the velocity distribution that is measured in 
the Itoh et al. A^-body simulations of clustering from an ini- 
tially Poisson distribution of identical particles that have no 
initial peculiar velocities. (These 10*-body simulations are 
described in more detail by Itoh et al. 1993). For ease of com- 
parison with published results, all velocities are in "natural 
units" (Saslaw et al. 1990; Itoh et al. 1993). The histogram 
in Fig. 4 shows the Press-Schechter model (eq. 25) when it 
is required to have the same dispersion as the simulations, at 
that epoch when the 10* particles have collapsed into Press- 
Schechter clumps, so that the spatial distribution of the par- 
ticles is described by equation (27) with b Si 0.75. (In the 
simulations, this occurs at an expansion factor a/ao = 8.) 
For comparison, the dotted line shows a Maxwellian that 
has the same dispersion as the histogram. The histogram is 
significantly different from the Maxwellian. 

The f{v) distribution that is measured in the simula- 
tions is represented by the dashed line, which shows the best 
fit of 

ft \^ a(T^(l-B) / 2/-, nx , 2n\at-2-l 

•^^"^ = r(a^^ + l) '-"^ {l-B) + av B) 

X g-a<x = (l-B)-a. = B ^^^^ 

to the simulation histogram. Equation (28) was first ob- 
tained by Saslaw et al. (1990), and has been shown to fit 
the velocity distributions in the Itoh et al. A^-body simu- 
lations with remarkable accuracy (Inagaki, Itoh fe Saslaw 
1992; Itoh et al. 1993). In equation (28), a and B are free 
parameters, and cr^ is the three dimensional mean square 



velocity. The top right corner of the plot shows the values 
of the parameters that provide the best-fit to the Press- 
Schechter histogram shown. These best-fit values are very 
similar to (within a few percent of) those measured in the 
A^-body simulations (Itoh et al. 1993). Similar agreement is 
obtained for all values of 6^0.8; at later times (6^0.8) it 
is not clear that these small (10*-body) simulations remain 
statistically homogeneous. Since the best fitting values of 
equation (28) to the A^-body f{v) are so similar to those of 
the Press-Schechter model (eq. 25), and in both cases the fits 
are very accurate, we conclude that equation (25) provides 
a good fit to the velocity distribution in the simulations. 

Since the Press-Schechter distribution evolves as the 
clustering develops, the model (eq. 25) can describe the evo- 
lution of the velocity distribution in the simulations. As the 
spatial distribution tends to the initial Poisson distribution 
(i.e., as 6 ^ 0), then equation (25) reduces to a Maxwellian, 
as required by linear theory. At later times, the distribution 
becomes significantly different from Maxwellian. This evo- 
lution simply reflects the changing fraction of particles in 
small clumps (small velocities for free particles that lose en- 
ergy to the expansion of the universe) relative to the fraction 
in massive virialized systems (having high virial velocities). 
Thus, as the clustering develops, the velocity distribution 
becomes skewed relative to the Maxwellian. This evolution 
is consistent with that measured in the simulations (Itoh et 
al. 1993). 

Given the accuracy with which equation (28) fits equa- 
tion (25), it is instructive to consider its functional form 
further. In effect, the parameter a simply rescales all ve- 
locities. Saslaw et al. (1990) hypothesize that the relation 
aV^ = Gm(l/r) = (Gm/a)(l/x), where m is the mass 
of each particle in their simulations, (l/x) is some mea- 
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Figure 5. Approximately self-similar evolution of the velocity 
distribution in A^-body simulations of clustering from an initially 
cold Poisson distribution. Curves show /(v/vrms) as a function of 
v/v-rms for a range of expansion factors. The bold solid line shows 
that equation (28), with B = 0.87, provides a good description of 
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sure of the average separation between particles measured 
in comoving coordinates, and a denotes the expansion factor 
since some fiducial time (usually taken to be the initial time, 
ao = 1), should always be satisfied. Since Gm{l/x) is some 
constant, their hypothesis means that acr^ = a'acr^, should 
be constant. This hypothesis appears to be in good agree- 
ment with what Itoh et al. (f993) measure in their A^-body 
simulations of clustering from cold Poisson initial conditions 
(see Model S in their Fig. ff; they use "natural units" in 
which their /? is proportional to our cr^ and their ^velocity 
is our B). Typically, the agreement with their hypothesis is 
good after the simulations are a few ( ^ 4) expansion factors 
larger than the initial size. 

Consider those epochs over which acr^ is constant. If _B 
is also constant at these epochs, then equation (28) evolves 
self-similarly. Fig. 11 of Itoh et al. (1993) shows that indeed, 
at late times, B is roughly independent of expansion factor. 
If we consider these Itoh et al. simulations (from Poisson 
initial conditions) as being roughly equivalent to simula- 
tions from Gaussian ra = initial conditions, then the Davis 
fe Peebles scaling solutions should apply to these (initially 
Poisson) simulations also. In this respect, the self-similar be- 
haviour (at times when acr^ and B are both approximately 
constant) of equation (28) can be thought of as approximat- 
ing the self-similar f{v) distribution that must obtain in the 
Davis fe Peebles description. 

Fig. 5 shows that, to a good approximation, the evolu- 
tion of /(?)) in these Itoh et al. simulations is self-similar. The 
histograms show f{v/vrnis) as a function ofv/v^ms when the 
radius of the simulation sphere is 4 (solid), 5.66 (dashed), 8 
(dot-dashed), 11.3 (dotted), 16 (dot-dot-dashed), and 22.63 
(solid) times the initial radius. The bold solid line that pro- 
vides a good fit to all the histograms shows equation (28) 



with B = 0.87; at each epoch it has been rescaled to describe 

So, equation (28) admits (at least approximately) self- 
similar behaviour, and this self-similar behaviour is in good 
agreement with the A^-body simulations. However, as shown 
in the previous section, self-similar behaviour arises nat- 
urally in the Press-Schechter model developed in this pa- 
per. Thus, the close agreement between the shapes of equa- 
tions (25) and (28) imply that the Press-Schechter model 
developed in this paper is able to explain the shape, and the 
self-similar evolution, of the distribution, f{v), of peculiar 
velocities. 

The close agreement between the distribution (eq. 28) 
obtained by Saslaw et al. (1990) and that considered here 
(eq. 25) is striking. However, there are two important dif- 
ferences between the Saslaw et al. (1990) derivation of f{v) 
and that presented here. Saslaw et al. used the assumption 
that fluctuations in kinetic energy are proportional to fluc- 
tuations in potential energy to derive equation (28) from the 
spatial distribution function (eq. 27). In particular, they as- 
sumed that (V'^) oc {GM/r), where the average is taken over 
randomly placed cells of some given size. They argued that 
this implies that {V^) oc {M). That is, they assumed that 
{Mir) was separable, that {Mir) = {Af){l/r), where the av- 
erage was taken over all space. However, our model of viri- 
alized, isothermal, Press-Schechter clumps, is constructed 
difl:"erently. In our model, {F^> oc {M^'^), where the average 
is taken over the Press-Schechter distribution of clump sizes. 
In this case, we showed that the averages are not separable. 

Secondly, since Saslaw et al. derive equation (28) for 
f{v) from equation (27) for f{N), their derivation suggests 
that b = B. In our model there is no reason for this equal- 
ity to hold, and, indeed, the simulations show that B > b 
always. In our model B is simply a free parameter that is de- 
termined by the best-fit of equation (28) to the A^-body or to 
the Press-Schechter f{v) histogram. Furthermore, whereas 
b increases as the simulation evolves. Fig. 5 suggests that, af- 
ter initial transients (due to the fact that in the simulations, 
particles have no initial velocities) have died out, then, to a 
very good approximation, B is constant, and the subsequent 
evolution of f{v) is self-similar. 

It is well known that equation (27) for f{N) is affected 
by random sampling (e.g. Lahav fe Saslaw 1992). On the 
other hand, equation (25) shows that f{v) in our model 
is the same for a subsample drawn randomly from some 
parent distribution as it is for the parent distribution. The 
simulations show that, indeed, f{v) is hardly changed by 
random dilution. In our model, therefore, the difference be- 
tween b [measured from f{N)] and B [from f{v)] can be 
used to provide information about the true underlying dis- 
tribution when only a randomly diluted subsample has been 
measured. 

Although Fig. 4 suggests that equations (25) and (28) 
are quite similar, they behave differently in the large v limit. 
These differences appear when f{v} ^ 10~* and so they are 
beyond the resolution of the 10*-particle simulations against 
which f{v) has been tested. Given the nature of the deriva- 
tion of equation (25), and its ad hoc treatment of the ve- 
locities between clumps in particular, we will not study the 
behaviour of this large v tail further. We simply note that 
the model presented here serves to illustrate that the distri- 
bution of velocities in the A^-body simulations can be un- 
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Figure 6. Example of the distribution function of pairwise relative peculiar velocities as a result of clustering from an initially Poisson 
distribution. Histogram shows this distribution when the expansion factor, a/ao, in the W-body simulations described in the text, is 8. 
At that epoch, the spatial distribution of all particles is well fit by equation (27) with b = 0.75. Dot-dashed line shows an exponential 
distribution that has the same dispersion as the histogram. Dotted line shows a Gaussian that has this same dispersion. Dashed line 
shows the distribution that is obtained by integrating equation (28) over directions perpendicular to the line of sight. Top right shows the 
value of the rms pairwise peculiar velocity, Ur rms, in "natural units" (Inagaki et al. 1992), as well as the values of the other parameters 
in equation (28), that are necessary to provide a good fit to the histogram. 



derstood by making rather simple assumptions about the 
properties and motions of Press-Schechter clumps. 

It is worth noting that the particular prescription we 
have chosen to specify the distribution of clump-clump ve- 
locities is not completely ad hoc. Firstly, it is clear that 
the clumps must have some net velocities with respect to 
each other. This is easily verified by locating the Press- 
Schechter clumps in the simulations, and then computing 
the net velocity of the center of mass of each clump. The 
shape of the distribution of clump velocities so obtained is 
better described by equation (28) than by, for example, a 
Maxwellian. However, we have shown that equation (28) is 
a good approximation to the distribution that is obtained 
by adding up a number of Maxwellian distributions with a 
(Press-Schechter) distribution of dispersions. In effect, our 
model assumes just such a prescription for clump-clump mo- 
tions. Since the motion of a clump of mass M is assumed 
to be drawn from a Maxwellian with dispersion related to 
M, on average in this model, more massive clumps move 
faster than less massive clumps. It is intriguing that mas- 
sive clumps have been measured to move faster than less 
massive clumps in simulations of a CDM (rather than an 
initially Poisson) universe (see caption to Fig. 5a in Zurek 
et al. 1994). 

Finally, recall that the distribution for f{v) that is de- 
fined by equation (28) is similar to that obtained from the 
Poisson Press-Schechter model developed in this section, for 
all measured values of b (as an example. Fig. 4 shows this to 
be true when b Si 0.75). However, notice that the prescrip- 
tion for obtaining f{v) from the Press-Schechter approach 
is very similar to that for obtaining f{u); both are integrals 
over Maxwellians having a Press-Schechter distribution of 
dispersions. The only difference between f{u) and f{v) is 



that, whereas the pairwise distributions are obtained by 
weighting the Press-Schechter mass functions by the number 
of pairs (oc M*''^ for truncated, singular isothermal spheres), 
the f{v) distribution in this section weights by the number 
of particles (oc M). If this difference in weighting terms is 
not significant, then, with the appropriate redefinitions of 
parameters, the distribution that is obtained by integrat- 
ing equation (28) over directions perpendicular to the line 
of sight should also provide a good approximation to the 
Poisson Press-Schechter distribution of pairwise velocities, 
/(«r; r). 

Fig. 6 shows an example of the line of sight pairwise ve- 
locity distribution that is measured in the Itoh et al. simula- 
tions. Since these simulations have only 10* (identical) par- 
ticles. Fig. 6 shows the distribution of line of sight, pairwise 
velocity differences for all pairs separated by less than some 
scale fcut . Although it includes contributions from pairs of 
varying separations, and so differs from the usual f{u;r) 
which is usually computed as a function of pair separation, 
it should provide a reasonable description of the shape of 
the more conventional f{u; r) on scales smaller than that 
of a typical clump. To insure that this is the case. Fig. 6 
is constructed with 2(3/47rrcut) = 300/9b, since clumps are 
assumed to have an overdensity of at least ITSph. 

The histogram shows this f{u; r) distribution when the 
simulation has expanded to a/ao = 8 times its initial ra- 
dius. That is. Fig. 6 shows the pairwise velocity distribution 
f{u; r) when the distribution of peculiar velocities them- 
selves, f{v), is given by Fig. 4. Dot-dashed and dotted lines 
show exponential and Gaussian distributions that have the 
same dispersion as the histogram. Whereas the exponen- 
tial distribution provides a much better fit to the histogram 
than the Gaussian, it has a higher, sharper peak than the 
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histogram. The reason for this sharper peak was discussed 
in the previous section. 

The dashed line that provides an even better fit to the 
histogram is obtained by integrating equation (28) over di- 
rections perpendicular to the line of sight. The top right of 
Fig. 6 shows the value of the rms pairwise peculiar veloc- 
ity, Mrrms, iu "natural units" (Inagaki et al. 1992), as well 
as the values of the other parameters in equation (28), that 
are necessary to provide a good fit to the histogram at this 
epoch. That the dashed line in Fig. 6 fits the histogram so 
well suggests that the Press-Schechter model of virialized 
isothermal clumps, when used to describe clustering from 
an initially Poisson distribution is, substantially, correct. 

Fig. 6 shows that the shape of the pairwise velocity dis- 
tribution for a distribution of particles that has evolved from 
an initially Poisson distribution, like that for a density field 
that was initially a scale-free ra = Gaussian random field, is 
well approximated by an exponential. Insofar as clustering 
from an initially Poisson distribution resembles clustering 
from an ra = Gaussian density field, this Poisson Press- 
Schechter model provides another test of the importance of 
discreteness effects. This is because an initially Poisson dis- 
tribution has a natural lowest mass scale — that of a single 
particle. Fig. 6 shows that discreteness effects may alter the 
small velocity part of the velocity distribution, whereas the 
large velocity tail is not strongly affected. These differences 
in the peak of the distribution are consistent with the discus- 
sion at the end of the previous section, regarding the effects 
of a small-mass cutoff in the mass distribution. 

Having established that the Poisson Press-Schechter 
model provides a good description of the pairwise veloc- 
ity distribution, and of the distribution of peculiar veloci- 
ties themselves, in (initially Poisson) A^-body simulations, 
we now consider the observations. Raychaudhury fe Saslaw 
(1995) show that the observed peculiar velocity distribu- 
tion function of a representative sample of galaxies (the 
Matthewson et al. 1991 spirals, and the Dressier et al. 1987 
ellipticals) is well described by integrating equation (28) over 
directions perpendicular to the line of sight. The close agree- 
ment between equations (28) and (25) (Fig. 4), and the fact 
that integrating a functional form like equation (25) over di- 
rections perpendicular to the line of sight gives a distribution 
that is nearly exponential (dashed and dot-dashed curves 
in Fig. 6), shows that the f{vr) distribution used by Ray- 
chaudhury fe Saslaw is well approximated by an exponential 
distribution. Thus, the Raychaudhury fe Saslaw results can 
be used to show that the exponential distribution provides a 
good approximation to the observed f{vr) distribution. This 
accuracy of the exponential distribution has also been noted 
by Bahcall, Gramann fe Cen (1994), in their study of the 
f{v) distribution of galaxies. It is intriguing that, as with 
the spatial distribution of galaxies, the approximately ex- 
ponential distribution of galaxy peculiar velocities is easily 
modelled by the Poisson Press-Schechter approach. 



4 DISCUSSION 

Section 2 used the Press-Schechter description of nonlin- 
ear clustering to obtain an expression for the distribution of 
pairwise relative peculiar velocities. The derivation assumed 
that Press-Schechter clumps are virialized, truncated, sin- 



gular isothermal spheres, that all have the same average 
density. Section 2 argued that if the velocity distribution 
within each Press-Schechter clump is Maxwellian, then the 
pairwise velocity distribution (when viewed along the line 
of sight and averaged over many different clumps) can be 
obtained by convolving Gaussians having a range of disper- 
sions. It showed that the distribution of dispersions can be 
related to the Press-Schechter distribution of clump masses. 
Fig.s 1 and 2 showed that this Press-Schechter model can 
explain the exponential shape of the line of sight distribution 
of pairwise velocities that is measured in A^-body simulations 
of gravitational clustering. This exponential shape also pro- 
vides a good description of the pairwise velocity distribution 
of galaxies. 

Section 2 also showed that, on small scales, the pairwise 
velocity dispersion in this model increases with scale. Since 
the size of a Press-Schechter clump, rvir, is only proportional 
to the cube root of the mass, the increase of the dispersion 
of f{u; r) with increasing pair separation, r, is expected to 
be small. The dispersions shown in Fig. 2 (corresponding to 
scales on which { ~ 10) are larger than those which obtain 
in the zero-separation limit (eq. 20). On these small scales, 
this scale dependence is in qualitative agreement with that 
suggested by the cosmic virial theorem (e.g. Peebles 1980; 
Efstathiou et al. 1988). 

The A^-body simulations show that as the separation 
between pairs increases and the pairwise velocity disper- 
sion increases, the distribution f{u; r) becomes slightly skew 
(Fig. 6 of Efstathiou et al. 1988). As the pair separation in- 
creases, the two members of any given pair are increasingly 
likely to be drawn from different clumps. These clumps need 
not be moving with the same net peculiar velocity, nor need 
they necessarily have comparable circular velocities. To date, 
there is no good model for clump-clump motions, so that the 
scale dependence of this skewness cannot be calculated ex- 
plicitly from the model presented in this paper. However, a 
skew distribution results when several Gaussians with dif- 
fering means and variances are added together. Thus, in 
our model, this skewness is easily explained. Moreover, it 
is possible to use the model presented here to constrain the 
clump-clump velocity distribution, though it has not been 
done here. This constraint follows from requiring that, in 
the limit of very large separations, when the two members 
of a given pair are almost certainly drawn from different 
clumps, then the line of sight pairwise velocity distribution 
must reduce to the form that is expected from models of 
bulk flows. 

Section 3 considered a simple model for the motions 
of Press-Schechter clumps. Using this model, it provided 
a description of the distribution of peculiar velocities f{v), 
rather than pairwise peculiar velocities f{u). Fig. 4 showed 
that the model provides a good description of the evolu- 
tion of the velocity distribution as an initially Poisson dis- 
tribution clusters gravitationally. Like the Press-Schechter 
model for the line of sight f{u), the line of sight f{v) dis- 
tribution is essentially a convolution of a Gaussian with the 
Press-Schechter distribution of dispersions, times a weight- 
ing term. The essential difference between the two distribu- 
tions is in this weighting term; it is oc Af (to account for 
the total number of particles) for f{v), and oc M"/^ (for the 
total number of pairs in an isothermal clump) for f{u). 

Given that the essential difference between f{u) and 
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f{v) is in the weighting term, the close agreement in Fig. 6 
between the dot-dashed curve (exponential distribution) and 
the dashed curve (integral of eq. 28 over directions perpen- 
dicular to the line of sight) has an interesting implication. 
Fig. 4 showed that equations (25) and (28) are extremely 
similar. Thus, the dashed curve in Fig. 6 is equivalent to 
integrating equation (25) over directions perpendicular to 
the line of sight. In other words, the dashed curve is a good 
approximation to the distribution that is obtained by sum- 
ming up Gaussian distributions with a Press-Schechter dis- 
tribution of dispersions, where each Gaussian is weighted by 
the number of particles rather than by the number of pairs. 
However, weighting by the number of pairs gives a distribu- 
tion that is similar to an exponential (Fig. 1; Section 2). So, 
the close agreement between the dashed and the exponen- 
tial (dot-dashed) curves in Fig. 6 shows that the shape of 
the integral over the Press-Schechter distribution of disper- 
sions is relatively insensitive to whether one is weighting by 
the number of pairs, or by the number of particles. 

The Appendix showed that, for the special case when 
clumps have a 'tophat' density profile, then weighting by 
the number of pairs is the same as weighting by the number 
of particles. For other density profiles, these two weighting 
functions will differ. Thus, Fig. 6 shows that the exponen- 
tial shape of the pairwise velocity distribution is primarily 
a result of the Press-Schechter distribution of dispersions, 
and is not so sensitive to the density profile within clumps. 
Alternatively, if the density profile is fixed, then one can 
group the weighting function with the Gaussian distribu- 
tion of (line of sight) velocities within clumps. In this case. 
Fig. 6 can be interpretted as evidence that the exponential 
shape will result even if the line of sight velocity distribu- 
tion within a given clump is not exactly Gaussian. In other 
words, the exponential shape of f{u) is not so sensitive to 
the exact density profile within clumps, nor is it particu- 
larly sensitive to the exact velocity distribution within each 
clump. Provided that the distribution of velocity dispersions 
within clumps is given by something like the Press-Schechter 
distribution, the exponential shape will result so long as the 
density profiles of clumps are something between tophats 
and isothermal spheres, and the velocity distribution within 
each clump is approximately Maxwellian. 

The distribution of dispersions is a consequence of as- 
suming that at any epoch, the circular velocity of a Press- 
Schechter clump may be related to its mass: Vc (xM^'\ This 
relation appears to be in good agreement with that measured 
in A^-body simulations (e.g. Bond fe Myers, in preparation). 
However, the Vc oc M^^^ relation is not consistent with the 
usual stability assumption that, after virializing, clumps, on 
average, do not change their physical size as the Universe 
continues to expand (e.g. Davis fe Peebles f977; Efstathiou 
et al. f988; Hamilton et al. f99f; Nityananda fe Padman- 
abhan f993). The reason for this inconsistency is clear. It 
arises from the fact that not all clumps of a given mass at a 
given epoch will have formed at that epoch; some may have 
formed at an earlier epoch. 

As a particular example, consider two clumps, each hav- 
ing the same mass, M, at a specified epoch, say, the present. 
Consider the case when one of these clumps was formed 
(i.e., virialized at f78 times the background density) at the 
present epoch, but the other formed at some earlier epoch 
(i.e., since the time when it virialized, it has survived with- 



out merging into a larger clump). By hypothesis, the clump 
that was formed earlier will have virialized at f 78 times the 
background density at the earlier epoch. However, the back- 
ground density decreases as the Universe expands. Since the 
background density at the time of formation determines the 
size, and so the circular velocity of the clump, then if the 
clustering is stable, the clump that was formed earlier will 
have a higher circular velocity than the clump of the same 
mass, M, that was formed at the present epoch. 

Thus, if the clustering is stable, then at any given epoch, 
there should be some scatter around the Vc oc M^/^ rela- 
tion. Since Vc oc p^^, it is clear that this scatter will not be 
large. Nevertheless, it may change the slope of the relation 
between the circular velocity and the mass. The discussion 
above suggests that as a result of this scatter, Vc oc M"^, 
where a^f/3, will be more accurate. Bond fe Myers (in 
preparation) find that a Si 0.29. Since this scatter is related 
to the probability that a clump of mass M at some earlier 
epoch is still of size M at some later epoch, it can be cal- 
culated self-consistently from the Press-Schechter approach. 
Thus, in principle, the Press-Schechter merger probabilities 
(Bond et al. f99f ; Lacey fe Cole f993) can be used to quan- 
tify the scatter around the Vc oc M^^^ relation. (In fact, the 
merger probabilities show that the probability of zero ac- 
cretion during some interval dt is zero. However, physically, 
one expects that a small amount of accretion in time dt will 
not change the density of the clump significantly. It is this 
notion, that a 'small' amount of accretion will not affect the 
density of a given clump significantly, which introduces the 
scatter around the density equals ITSph relation.) 

Since this scatter is small (recall that Vc oc p\J^), while 
it may affect the scale dependence of the distribution of pair- 
wise relative velocities, it is unlikely to affect the shape of 
this distribution. Therefore, despite the minor inconsistency 
between the stable clustering prescription and the Press- 
Schechter model for the circular velocities of clumps, the 
conclusion that the Press-Schechter model can easily ex- 
plain the exponential shape of the distribution of pairwise 
relative velocities along the line of sight remains valid. 

There is another minor discrepancy between the stable 
clustering hypothesis and the Press-Schechter model dis- 
cussed above. It arises from the fact that we assumed that 
all Press-Schechter clumps were isothermal spheres, inde- 
pendent of the slope n of the initial power spectrum. In fact, 
stable clustering and the scaling solutions of the BBGKY hi- 
erarchy (Davis fe Peebles f977; Pebbles f980) require that, 
in the highly nonlinear regime, the density within a Press- 
Schechter clump should scale as r~'\ with 7 a function of 
ra, rather than as as assumed in Section 2. However, as 
we have argued in this section, the exponential shape of the 
pairwise velocity distribution is relatively insensitive to the 
density profile within clumps. Thus, a more exact calcula- 
tion that incorporates the n dependence of the (small scale) 
density profiles within clumps will not change our conclu- 
sion that the Press-Schechter model can easily explain the 
exponential shape of f{u). 

Although this paper has focussed on self-similar models 
for the initial density field (and so on cosmological models 
that have critical density), equation (5) is quite general; it 
may be used with any Press-Schechter multiplicity function. 
So, in principle, it may be used to compute the pairwise rel- 
ative velocity distribution for cosmological models in which 
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the density parameter is less than unity. As reliable pecu- 
liar velocity measurements on small nonlinear scales become 
available, equation (5) may provide a way to estimate Q. 

However, there is a complication when comparing the 
Press-Schechter models of clustering from arbitrary initially 
Gaussian density fields with the observations. It arises from 
two uncertainties. The first is simply that observations mea- 
sure velocities of discrete particles (galaxies), and it is not 
clear how to relate the Press-Schechter distribution of pair- 
wise velocities in a continuous density field to that measured 
for discrete particles. For this reason Section 2 used a Pois- 
son sampling argument to show that the effects of discrete- 
ness could be calculated self-consistently for the Efstathiou 
et al. A^-body simulations (Fig. 3). In addition. Section 3 
studied the Press-Schechter description of clustering from 
an initially Poisson distribution — since an initially Poisson 
distribution has a natural lowest mass scale. It argued that, 
insofar as clustering from an initially Poisson distribution 
should resemble that from an initially Gaussian, scale-free, 
ra = model, the Poisson Press-Schechter model can be 
treated as a self-consistent way of imposing the small-mass 
cutoff. The similarity of the Poisson Press-Schechter model 
to the ra = Gaussian model, and the accuracy of the Pois- 
son sampling model for models with other values of ra sug- 
gests that, at least for the particles in the simulations, the 
effects of discreteness can be understood. Except for the 
peaks of the distributions, the overall shapes of the pairwise 
velocity distributions are hardly effected by the discreteness. 

The second difficulty in relating these Press-Schechter 
models to the pairwise velocity distribution of galaxies is 
the well known problem of biasing. Whereas the simulations 
trace the evolution of dark matter, it is not clear that the 
baryonic matter from which (luminous) galaxies are thought 
to form, necessarily has the same distribution as the dark 
matter. In particular, whereas the A^-body simulations sug- 
gest that the Press-Schechter model of virialized isother- 
mal spheres is an adequate description of the distribution 
of sizes of dark matter clumps, it is not clear that it can 
also describe the distribution of sizes of galaxy clusters, or 
the distribution of galaxies within a cluster. If the number 
of galaxies associated with a given dark matter clump is re- 
lated to the mass of that clump, and if the distribution of 
galaxies within a cluster traces the distribution of the dark 
matter, then the Poisson sampling model provides a way to 
predict the velocity distribution of galaxies from the Press- 
Schechter description of the dark matter. If the biasing is 
more complicated, then it is not clear that the shape of the 
velocity distribution of the dark matter need have any simi- 
larity to that of the galaxies. Therefore, it is intriguing that 
the line of sight pairwise velocity distributions of the dark 
matter distributions studied in this paper are, like the ob- 
servations presently available, very well fit by exponential 
distributions. 
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APPENDIX A: THE DISTRIBUTION OF PAIRS 

This Appendix computes the number of pairs as a function of scale, s, within a sphere of mass M and radius R. Let N{s) ds 
denote the number of pairs within a sphere of radius R and mass M that have separations in the range ds about s. Then 
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Of course, N{s) = for all s > 2R. Clearly, N{s) depends on the density distribution within the sphere. For a sphere of 
uniform density, p{x} = 3M/4:wR^ is independent of position, so 
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The total number of pairs in the clump is 
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as expected. The factor of 1/2 is included to avoid counting pairs twice, and the upper limit on the integration is set by the 
fact that N{s) = for separations larger than the diameter of the clump. The useful result is that, to lowest order in s/R, 
the number of pairs scales as M'^ / R^ . Notice that, for a uniform sphere, M R^ . This means that the number of pairs scales 
as the number of particles M. 

The calculation for a truncated singular isothermal sphere of mass M and radius R can be performed analogously. The 
density, p, within a singular isothermal sphere scales as where r is the distance from the center of the sphere. So, let 

p = A/r'^ . The constant A is set by requiring that the integral of the density over the sphere equal the mass, M. Then 

fR j"K l'2'K l'2'K t'^'^ 

dr sm9d9 d<t> I I 
lo Jo Jo Jo Jo 

2 {'R—s {'TT {'2tt {'2tt 

drsm9d9d4>\ / ^r—. — + 



N{s) oc 



M ^ ^ 
47r_R 



M 
AttR 



M 
AttR 

AttR. 



Jo 
2 i-R 



s'^ sin9' d9' d(l)' 
^2 _|_ ^2 _ 2sr cos 9' 

" sin 61' d6i' d,/.' 



Jo 



_|_ j,2 
2n /-e' 



dr sin 9 d9 dct> 



2sr cos 9' 

s^ sm9' d9' d<f>' 



g2 _|_ ^2 



2sr cos 9' 



R- 



dr , ( s -\- r 
— In 

r \s — r 



R 



R 



dr , 
— In 

r \s — r 



(A6) 



This must be evaluated for two separate cases. When s > R — s, oy equivalently, when R > s > R/2, then 
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When s < _R — s, or, equivalently, when s < -R/2, then 
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In either case, to lowest order, N(s) oc M'^ /R^ . 

The text considers Press-Schechter clumps, which are spheres that all have the same average density, so that M (x . 
This means that, to lowest order, the number of pairs scales as M'^/R^ oc M for a uniform sphere, and as M^/R^ oc M"/^ for 
truncated singular isothermal spheres. This is the result used in the text. 

In this paper we are primarily interested in the way in which the number of pairs of a given separation, s, depends on 
the mass, M, and radius, R, of a given spherical clump. However, the calculations in this Appendix are also useful for the 
following problems. Suppose we are interested in calculating the two-point correlation function for particles within a sphere. 
To do so, we need to be able to estimate the number of pairs as a function of scale, under the assumption that the particles 
are distributed uniformly at random throughout the sphere. If the sphere were of infinite radius, and if the density within the 
sphere was p = 3M/4:wR^ , then the number of particles that are a distance s from a given particle is proportional to Awps'^, 
so that the total number of pairs that have separation s is half of pV x Awps'^ = (3M/R^) . However, if the sphere were of 
finite size, then the number of pairs as a function of scale is modified. The correct expression for N(s) for the uniform sphere 
case is equation (A4). In the limit as R ^ oo, it reduces correctly to the infinite sphere case. For separations that are on the 
order of 10% of the size of the sphere, the correction to the infinite sphere case is on the order of 10%. 

Next, suppose that we are at the center of a sphere of approximately uniform density, and again, we are interested in 
estimating the correlation function for the particles within the sphere. Suppose, however, that our efficiency at detecting 
partices, our selection function, is a function of the distance to the particle. For convenience, suppose that that the efficiency 
decreases as 1/r^ . Then the calculation in this Appendix provides a way of estimating the number of pairs as a function of 
separation, assuming that the underlying distribution is uniformly random. In this case, p{r) plays the role of the selection 
function. 



